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Cn ' Abstract. In this paper we investigate a variant of the classical mixmaster 

» I ' universe model of anisotropic cosmology, where the spatial sections arc non- 

C^ commutative 3-tori. We consider ways in which the discrete dynamical system 

describing the mixmaster dynamics can be extended to act on the noncommu- 
tative torus moduli, and how the resulting dynamics differs from the classical 
f^ ' one, for example, in the appearance of exotic smooth structures. We discuss 

properties of the spectral action, focussing on how the slow-roll inflation po- 
tential determined by the spectral action affects the mixmaster dynamics. We 
relate the model to other recent results on spectral action computation and 
O ' we identify other physical contexts in which this model may be relevant. 
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1. Introduction 

Noncommutative geometry has been used extensively as a method for the con- 
struction of models of particle physics, cosmology, and gravity coupled to matter. 
The point of view adopted in these models is conceptually similar to the extra di- 
mensions of string theories, in the sense that one replaces a 4-dimensional spacetime 
manifold with a product (or fibration), where the fibers are the "extra dimensions". 
However, unlike in string theory, these extra dimensions are not manifolds but non- 
commutative spaces. Moreover, on this product geometry encompassing spacetime 
directions and noncommutative extra dimensions, one has an action functional, the 
spectral action^ which is a natural action functional for pure gravity in noncommu- 
tative geometry. The key idea is that pure gravity on a noncommutative space that 
is the product of a spacetime manifold and a suitable noncommutative fiber looks 
like gravity coupled to matter from the spacetime point of view. More precisely, the 
asymptotic expansion for the spectral action functional delivers terms that contain 
gravity terms such as the Einstein-Hilbert action and a Lagrangian for the matter 
content of the model, which depends on the choice of the finite geometry. For a 
simple choice of a finite dimensional algebra as the fiber noncommutative space one 
recovers the Standard Model Lagrangian [TJ , extensions of the Minimal Standard 
Model that include right handed neutrinos |12| and supersymmctric extensions [8]. 

Unlike other high-energy physics models involving noncommutativity, such as 
those arising in some string theory compactifications, there is here no noncom- 
mutativity in the spacetime coordinates, only in the extra dimensions. Moreover, 
these particle physics models live naturally at unification energy. While there is, 
at this point, no clear picture of how the models should be extended to higher 
energies, it has been frequently proposed that, when moving towards the Planck 
scale, the model should incorporate more noncommutativity, which will eventually 
involve the spacetime coordinates as well. Moreover, one expects that more "se- 
riously noncommutative" spaces (that is, not Morita equivalent to commutative 
ones) should appear as one approaches the Planck scale. In this paper we do not 
attempt to answer the question of how to extend the noncommutative geometry 
models of gravity coupled to matter towards the Planck scale, but we describe a 
toy model case for what a behavior of the type suggested would look like, in a 
geometry that is at the same time sufficiently simple to be explicitly computable, 
but sufficiently complex to exhibit a nontrivial behavior. 

Our model is constructed by adapting a very well known example of classical 
cosmologies that exhibit a chaotic behavior, namely the mixmaster cosmological 
models of [5] , [26] , based on the Kasner metrics and on a discrete dynamical system 
related to the continuous fraction expansion, that governs the succession of mix- 
master cycles and Kasner epochs. We input the noncommutativity in this model 
by replacing the spatial sections of this cosmology by noncommutative tori. 



NONCOMMUTATIVE MIXMASTER COSMOLOGIES 3 

In ij2] we recall the classical mixmaster universe model, which we formulate in 
the case where the spatial sections are 3-dimensional tori T^. We recall the main 
properties of the discrete dynamical system that models the mixmaster dynamics, 
and its relation to the Kasner metrics and the Kasner epochs and cycles of this 
anisotropic and chaotic universe model. 

In iJ3]we recall some basic properties of noncommutative 3-tori, as noncommuta- 
tive algebras and as spectral triples (noncommutative Riemannian manifolds). We 
describe two different possible ways to extend the mixmaster dynamical system to 
act on the moduli of the noncommutative tori and not only on their metric struc- 
ture. We show that one of these models leads naturally to the occurrence of exotic 
smooth structures in this noncommutative cosmological model. 

In 21 we discuss inflation scenarios derived from the spectral action functional 
and we construct a toy model of deformations of parameters that gives rise to 
a transition from an early universe noncommutative mixmaster cosmology to a 
commutative inflationary cosmology. 

Finally, in ^ we discuss other aspects of the model, related to properties of the 
spectral action, and we identify physical settings in which this type of model may 
be relevant. 

2. The mixmaster universe 

We review briefly in this section the basic properties of the mixmaster universe 
models in general relativity, and some aspects of the geometry that wc need for the 
noncommutative generalization that follows. 

2.1. Kasner metrics. The mixmaster universe models were developed (see [5], 
[26]) as interesting cosmological models exhibiting strong anisotropy and chaotic 
behavior. One considers anisotropic metrics 

(2.1) ds^ = -dt^ + a(tf dx^ + b{tf dy^ + c{tfdz^, 
with the functions a{t), b{t), c{t) of the Kasner form 

(2.2) ds^ = -dt^ + t^P'dx^ + t^P^dy^ + t^P^dz^, 

where the exponents Pi,P2,P3 satisfy the constraints J^iPi = 1 ^^^ J2iPi ^ 1- 

2.2. Discrete dynamical system. The mixmaster universe is a solution of the 
Einstein equation, built out of approximate solutions that look like Kasner met- 
rics for certain intervals of time (Kasner eras), with a discrete dynamical system 
governing the transition from one era to the next in terms of the change of Kasner 
exponents in the metric, see [531 ^^'^ E3' ESI- One sets 

—u 
Pi = 

(2.3) P2 = 

P3 = 



1 + 1 
1 



1 + u + u^ 
u{l + u) 



1 + u + u^ 

The dynamics is discretized by starting, at the beginning of each Kasner era, with 
a value u„ > 1 of the parameter u. Each era is then subdivided into shorter cycles, 
determined by decreasing values u„, m„ — 1, u„ — 2, . . .. Within each of these cycles 
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the metric is approximated with a Kasner metric p.2|) with exponents (|2.3|) with 
fixed u = Un — k. This sequence of Kasner cycles stops when the next value Un — k 
becomes smaller than one (but still positive). Then one passes to the next Kasner 
era with the transformation u i— > 1/w and restarts the sequence of Kasner cycles 
from this new value. Thus, the transformation of the parameter u that marks the 
passage from the beginning of one Kasner era to the the beginning of the next is 
the well known dynamical system 

(2.4) T -.Un^ Un+l = F T, 

Un - [Un\ 

which is the shift of the continued fraction expansion Tx = l/x— [l/x] with x„_|-i = 
Txn and u„ = l/a;„. 

Moreover, at the start of each new Kasner era, a permutation of the three spatial 
directions occurs, which reassigns the role of the direction responsible for the main 
dilation or contraction and of the two oscillating directions. As shown in [27| . 
[25] , in terms of the dynamical system (|2.4[) and the shift of the continued fraction 
expansion, this permutation can be described in the following way. Identify the 
three spatial directions with the three points of P^(F2) via 1 M- a;; oo n- y; i-^ z. 
Then the shift Tx ~ l/x— [l/x] of the continued fraction expansion on [0, 1] extends 
to a discrete dynamical system on [0, 1] x P-'^(F2) by 



(2.5) T:{xr ' ^ 



-[l/x] 1 
1 



This dynamical system (and its invertible two-sided extension) is in fact the dy- 
namical system that gives the coding of geodesies on the modular curve Xq{2) = 
EI/ro(2), as shown in [57], [55], which leads to the description of solutions of the 
mixmaster universe dynamics in terms of geodesies on this modular curve, as re- 
called in Section [2731 below. 

Clearly, the description of mixmaster universe cosmologies in terms of this dis- 
crete dynamical system only leads to the construction of an approximate solution of 
the Einstein equation, and one can then argue with more subtle analytic methods 
in what sense this approximate solution is close to an actual solution, but we will 
not be dealing with the approximation problem in this paper. 

2.3. Mixmaster data and geodesies on the modular curve Xq{2). It was 
shown in |27] , |28] , that the solutions of the discretized mixmaster dynamics are in 
one-to-one correspondence with geodesies on the modular curve Xo(2) = IHI/ro(2), 
with To (2) C SL2(Z) the congruence subgroup of level two. 

Recall from [27], [28] that every infinite geodesic on Xq{2) not ending at the 
cusp can be coded by data (a;,s) = (w~,w+,s), with uj"^ G P^(M) and s £ P^(F2), 
and where u!^ can be chosen with a;+ G [0, 1] n (M \ Q) and uj~ G (— c)o, — 1] H 
(K \ Q). These can be written in terms of continued fraction expansion as a;+ = 
[fco, fci, fc2, . . . , fc„, . . .] and uj^ = [/c-i; fc-2, • • ■ 7 fc-n, • • ■] and are acted upon by the 
shift as 
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Geodesies on Xo(2) are parameterized by the orbits of the data (w, s) under the 
action of the shift T. 

The data (w, s) in turn determine a solution of the mixmaster dynamics, by 
assigning the \un] to be the digits fc„ of the continued fraction expansion of a; and 
the alternation of the spatial directions being determined by the action of T on 
the element s G P^(F2), according to the identification mentioned above between 
points of P^(F2) and spatial axes. 

2.4. Mixmaster tori. Observe that the metric ()2.2|) can be considered equally 
on a spacetime whose spatial sections are topologically a flat space M'^ or whose 
sections that are topologically tori T^ . Wc will focus on the latter possibility 
for our noncommutative model. Thus, wc assume that the spacetime manifold is 
topologically a cylinder T^ x R, endowed with the Lorentzian metric of the Kasner 
form ([221). 

In the case of a torus T'^, the evolution in time {Tf,gt) with the Kasner metric 
g^ = t^P^dx^ + t^P^dy^ + t^P^dz^ has volume Vo1(T3) = tYol{T^), since Yol{T^) = 
^pi+P2+p3Yq^2^3j^ where we are assuming that pi + P2 + Ps = 1- 

Consider the Dirac operator ^t on {T^,gt), associated to a choice of a spin 
structure s on T^. On the 3-torus T^ there are eight different spin structures Sj. 
We recall the following result (see [2], [23]) on the Dirac spectrum. 

Let T^ ~ R-^/A be a 3-dimensional torus, with A a lattice in R'^. Let {ti, T2,ts} 
be a basis for A and let A^ be the dual lattice with dual basis {'Wi,W2,W3}. The 
eight spin structures arc classifies by eight vectors {s = (51,52,53) |Si G {0,1}}, 
where the value of each 5; distinguishes whether the spin structure on each of the 
three directions Vi is twisted or untwisted. In fact, on the circle S^, the spinors for 
the two possible spin structures can be identified with 

T{S\S) = {i/^ : R ^ C I iit + 27r) = ±iit)}, 

and the Dirac operator —i-i has cigcnfunctions 4'k{t) — exp(2fc7rit) for the trivial 
spin structure and ipkit) = exp((2/c + l)Trit) for the other one. 
On T^ — R'^/Z'^, the Dirac operator ^ is of the form 

(2.6) ^ = -i{<Ji— + cr2j-+a3 — )^ 

ox oy oz 

More generally, on T^ = M'^/A, one can write the Dirac operator in the form (up 
to a possible overall additive shift) 
3 

(2.7) ^ = -z ^(r, • d)<jj = -z((ti • d)ai + (t2 • d)a2 + (rg • d)a3) 

where d = {81,82, 83) = (^, ^, ^), 8rj = tj ■ 8, and aj are the Pauli matrices. 

Then on the 3-torus T^ = M'^/A, the Dirac operator ^ on the spin structure 5 
has spectrum 

1 ^ 

(2.8) Spec(^) = {±2tt\\w+-'^SjWj\\ |u;e A''}. 



»■ ^ 


--9--iA. 
dy dx 


~*al 


d ,■ d 


''dz 
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In particular, in the case of the torus Tj^ = M.'^/At, where the dual lattice A^ is 
spanned by the basis of vectors {t'''^ei,tP^e2,t'''^e3}, with ei the standard orthonor- 
mal basis, the Dirac operator is given by 

(2.9) ^--(-*-^^|: + -^^-^^| + -3t-|) 

/ -it-p^4- -t-p^ ^ - it-p^ ^ \ 

1 oz oy ox \ 

^\ t-p^- 4- - it-p^ 4- it-p^-S- ' 

\ ay ox oz I 



The spectrum of ^t, for the spin structure s, will then be of the forin 

1 ^ 
(2.10) Spec(^t) = {±27r||(i-Pifc,rP^m,i-P3n) + - ^Sji-P^ejH | (fc,m,n) G 1?}. 



2.5. Mixmaster dynamics on classical tori. The dependence of the Kasner 
exponents on the w-parameter as in (|2.3p and the discrete dynamical system (j2.4p . 
with the permutations (j2.5p . correspondingly determine a sequence of 3-tori T^ 
with Dirac operators as in (|2.10p . where the exponents pi are functions of m = w„ 
through p.3p and the permutation ()2.5p of the coordinate axes. 

As in Section 12. 3[ passing of each Kasner era is marked by the transition w„ t-> 
w„+i and Sn ^^ Sn+i, where w„ = Xjxn and x„+i = Ta;„, and with s„ g P-'^(F2), 
with s„+i = Tsn with the map T as in p.Sp . In particular, the action of T that 
gives the transition from the nth to the (n + l)st Kasner era is given by the action 
of the matrix 

(2.11) 7n = ( "f""^ \ ) e GL2(Z), 

which simultaneously acts on changing the metric of the torus and on permuting 
its three generators. 



3. NONCOMMUTATIVE 3-TORI AND MIXMASTER COSMOLOGIES 

We now reconsider the model of mixmaster universe described above, in a setting 
where the three spatial coordinates that give the 3-torus T^ are replaced by a non- 
commutative 3-torus Tgi . We describe these noncommutative spaces as Riemannian 
geometries, in the noncommutative sense, that is, as spectral triples. We then pro- 
pose different possible ways in which the discrete dynamical system that describes 
the evolution of a classical mixmaster cosmology can be extended to involve also 
an action on the parameters (the moduli) of the noncommutative torus itself, so 
that not only the metric structure, but also the underlying noncommutative space, 
evolves along with the succession of Kasner epochs. 

3.1. 3-tori as spectral geometries. In noncommutative geometry the analog of 
a Riemannian spin manifold is described by the data of a spectral triple. These 
consist of a triple {A,'H,D), where A is an involutive, dense subalgebra of a C*- 
algebra closed under holomorphic functional calculus, together with a (faithful) 
representation tt : A —> B{T-L) as bounded operators on a separable Hilbert space 
%, and a "Dirac operator" D. The latter is a self-adjoint (unbounded) operator, 
densely defined on K, with compact resolvent and satisfying the condition of having 
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bounded commutators with the elements of the algebra, [TT{a),D] G B{H), for all 
a e A. 

A smooth compact Ricmannian spin manifold X is a special case of a spectral 
triple, where the data (A,n,D) are given by ^ = C°°{X), U = L'^{X,E), with S 
the spinor bundle, and D = ^x the Dirac operator. 

Thus, the mixmastcr tori T^ — M.^/At described above are spectral triples with 
A = C°°{T^), and with H = L'^{T^,§s), where Sg is the spinor bundle for the spin 
structure s ~ (51,52,53) as above, with 

(3.1) L\T^,S,) = {i':R^ ^C^li'eL^, ip{x + tP^Vj) = {-iy^iP{x)}, 
for X ~ {x, y, z) G M.^ and {tj} the basis of A, and with Dirac operator 

3 

(3.2) ^ = -2Y,t-P'irrd)<yj = -^it-'"iri-d)a,+t-P-iT2-d)a2 + t-P-iT3-d)a3). 

Here we scaled the spatial coordinates Xi 1— > t^'Xi, so that dxi H> t^^^dxi, ^-'^d fixed 
the Ti, instead of scaling A n> A^ as before: the resulting p.2p is the same. 

3.2. Noncommutative 3-tori and their spectral geometry. A noncommuta- 
tive 3-torus is the universal C*-algebra ^e generated by three unitarics Ui, U2, U3 
with the relations 

(3.3) U,Uk = exp{2nie,k)UkUj , 
where Q = (Ojk) is a skew-symmetric matrix 

/ 03 -62 

(3.4) e = -03 01 

V 02 -01 

It can also be described as the twisted group C* -algebra C*{1? ^gq), where the 
f/(l)-valued 2-cocycle tre is given by (see [H]) 

(3.5) cre{x,y) = exp{TTi{x,&y)). 

A detailed discussion of the main properties of noncommutative 3-tori can be found 
in a. 

Recall that, for a noncommutative 3-torus, the algebra of smooth functions is 
defined as 

(3.6) A^^{XeAe,X^ ^ a™,„,fe C/rC/^'C/3'^ | a = (a™,„,fc) G 5(Z3,C)}, 

m,n,fctzZ 

which is the noncommutative analog of functions on the 3-torus with rapidly de- 
caying Fourier coefficients. 

It was recently proved by Venselaar in |41j that all equivariant real spectral 
triples on rank n noncommutative tori are isospectral deformations, in the sense of 
|15j . where spin structures are Dirac operators on commutative flat n-dimensional 
tori T". The case of rank- two tori was previously shown in [37j. 

Thus, these spectral triples will all be of the form (Aq , L^{T^, §), ^), where § is 
the spinor bundle for one of the eight spin structures on the ordinary torus T^ and 
the Dirac operator ^ is of the form p.7p (up to an overall additive constant) with 
{r,}j=i,2,3 a basis. 
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3.3. Noncommutative 3-tori, moduli, and mixmaster evolution. Now we 

consider again tlie discrete dynamical system of mixmaster evolution, where the 
passage from one Kasner era to the next is determined by the action of the matrix 
7n G GL2(Z) of (j2.1ip . In this noncommutative setting, in addition to the action 
on the exponents of the Kasner metric and the permutation of the spatial directions 
(here given by a permutation of the three generators of the noncommutative torus 
algebra) and by the corresponding action on the Dirac operator in (|3.2|) as in 
the commutative case, one also has the modulus Q of the noncommutative torus, 
which is trivial in the commutative case. Thus, one can propose extensions of the 
mixmaster dynamics where the modulus Q is also acted upon in the transition from 
one Kasner epoch to the other. We present in the following subsections examples of 
two possible such extensions of the mixmaster dynamics and we illustrate some of 
their properties. As we discuss below, these will have some interesting consequences 
on the properties of the noncommutative mixmaster cosmologies that differ from 
their classical counterparts. 

3.4. Moduli evolution by auxiliary mixmaster data. In this scenario, we 
assume given a choice of mixmaster data, by which we mean a choice of a geodesic 
on the modular curve Xo{2), or equivalently a choice of data (w, s) up to the action 
of the shift T as recalled in Section 12.31 This means that we have an assigned 
sequence of matrices 7„ G GL2(Z), of the form (|2.11|) . We also have a sequence s„ 
of elements in P^(F2), which determines, in each Kasner era, which spatial direction 
is the dominant direction driving expansion or contraction. 

Recall then that, if we write the modulus for the noncommutative 3-torus as 
the vector 9 = {9i, 62,03) G K"^ of the three parameters out of which the skew- 
symmetric matrix 8 of p.4p is built, then for any matrix M € GL3(Z) we have an 
action 9^ ~ M0 and we let Q' be the skew-symmetric matrix corresponding to the 
new values 6 . 

The fact that, in each era of the era of the mixmaster dynamics, there is one of 
the three spatial directions that dominates the contraction or expansion, identified 
by the given sequence s„ allows us then to define an action of the matrices 7„ G 
GL2(Z) on the parameters of the noncommutative torus, by embedding GL2(Z) 
inside GL3(Z), so that it acts as the identity on the parameter associated to the 
dominant direction and as 7„ on the other two. Since this action is accompanied by 
a permutation of the directions, at the following change of Kasner era the dominant 
direction will have changed and the embedding of GL2(Z) inside GL3(Z) used in 
defining the action on 9 will change accordingly. 

Thus, for example, if at the nth Kasner era the first coordinate is the dominant 
direction, we obtain the transformation 

(3.7) 9 = {91,62,93) ^ 9! = {9,, -[u^^]92 + 93,62). 

As we discuss in the next subsection, the choice of this model for mixmaster evo- 
lution of the noncommutative torus moduli has an interesting consequence: along 
the sequence of Kasner epochs, not only the metric structure of the noncommutative 
tori undergoes a sequence of transformations analogous to the classical mixmaster 
dynamics, but at the same time the smooth structure of the noncommutative 3-torus 
undergoes a sequence of transitions to different exotic structures. 
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3.5. Smooth structures and Kasner eras. In ordinary commutative geometry, 
the first occurrence of exotic smootli structures, meaning examples of smooth man- 
ifolds that are homeomorphic but not diffcomorphic, occurs in dimension four. In 
noncommutativc geometry, however, the simplest example of exotic smooth struc- 
tures is known to occur already in three dimensions, for noncommutativc 3-tori. 

Wc consider here the algebra A^ of smooth functions described in p.6p . 

One can then introduce the following two equivalence relations on the modulus 
e (see 13): 

• e - 9' ^ ^' = Aie, with M e SL3(Z); 

• e « 9' ^ ^' = M0, with M G GL3(Z). 

Then (see j^) one has algebra isomorphisms: 

(3.8) ^0-^0' ^ 9 « 9' while A^ ~ A^, ^ 9 - 9'. 

The C*-algebra case follows from [6], [4] and [21], while the smooth subalgebras 
case uses [6], [7], [16]. 

One sees from this result that the noncommutativc 3-tori have exotic smooth 
structures: any two tori with 9 w 9' through a matrix 9^ = M0 with M S GL3(Z) 
but M ^ SL3 (Z) are homeomorphic (in the sense that the algebras of continuous 
functions are isomorphic) but not diffcomorphic (in the sense that the algebras of 
smooth functions are not isomorphic). 

Thus, an interesting phenomenon happens in the mixmaster dynamics, whereby 
the passage to each successive Kasner era, which is determined by the action of a 
matrix M G GL3(Z), which has det(A/) = det(7„) = —1, changes the torus Tg, to 
a new torus Tq, which is homeomorphic, but with a different smooth structure. 

The topic of exotic smooth structures and their relevance to physics has been 
explored in various aspects in recent years, see for instance [1], [17]. This simple 
observation about the noncommutativc tori shows that, when allowing noncom- 
mutativity in the space coordinates, one can more easily encounter phenomena 
involving exotic smoothness, such as, in this case, changes of smooth structure. 

3.6. Moduli evolution by internal mixmaster data. We propose here another 
possible way to extend to the torus moduli the discrete dynamical system defining 
the mixmaster dynamics. In this case, one does not assume a given classical mix- 
master solution, but constructs it directly in terms of the torus moduli themselves. 

In this case, to define the action of the mixmaster dynamics on the modulus 9, 
we recall the following equivalent description of the noncommutativc 3-torus ^ej 
see [4]. One can view the 3-dimcnsional noncommutativc torus ^e a-s a crossed 
product C*-algebra for an action of Z on a 2-dimensional noncommutative torus: 

(3.9) ^0=A3>^aZ, 

where Ae^ is the 2-dimcnsional noncommutativc torus generated by two unitarics 
U and V with the relation VU = e^''''^''UV, and the action a : Z ^ Aut(^03) is 
given by 

a{U) = e^^'^'^^U, a{V) = e-'^'"''^V. 

One can then change the parameters of the noncommutative 3-torus, in passing 
to the next Kasner era, by acting on the 2-dimensional noncommutative 2-torus 
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Ae^ by a Morita equivalence, implementing the change of parameter given by the 
action of the matrix 7„ on ^3 by fractional linear transformations 

(3.10) e, ^ j,,{e,) = ~[""f^ + ^ = 1 - M, 

for 7„ as in (j2.1ip . but where now the integers A:„ = [u„] are the digits of the 
continued fraction expansion of 63 itself. 

In the particular case where the parameter 63 is a quadratic irrationality, namely 
an irrational number that is contained in some real quadratic field embedded in R, 
then the digits of the continued fraction expansion of 63 are eventually periodic, 
and there is a natural choice of the mixmaster data {uj,s), with u> = {^3,^3}, 
with 6*3 the Galois conjugate of 6^. This choice corresponds to a closed geodesic in 
^o(2). 

Notice that here we are acting by Morita equivalences of the 2-dimensional non- 
commutative torus Ae^ , which are implemented by an action of GL2 (Z) , while for 
3-dimensional smooth noncommutative tori in the generic case, the Morita equiv- 
alences arc implemented by the action of the group S'0(3, 3|Z), see [12], [33] and 
[20] for a complete classification up to Morita equivalences. 



4. The spectral action and inflation scenarios 

We now consider the action functional, the spectral action, for the noncommu- 
tative mixmaster cosmologies. 

The spectral action functional is defined as Tr(/(£'/A)), where A is the energy 
scale and / is a test function, usually a smooth approximation of a cutoff function. 
This is regarded as a spectral formulation of gravity in noncommutative geometry. 
This action functional has an asymptotic expansion at high energies (see (|5.10p 
below). Thus, one can approach the computation of this action functional either 
by explicit information on the spectrum and a computation of the series defined by 
Tr(/(Z3/A)) (as in [TU], [H], [H], [SI], [32), or else through its asymptotic expansion 
and the computation via heat kernel methods and local expressions in curvature 
tensors of the various terms in the expansion, as, for instance, in [9], [12], [22], [24] . 

The computation of the terms in the asymptotic expansion shows (see for in- 
stance [12]) that one recovers the terms in the usual classical action for gravity, 
namely the Einstein-Hilbert action and the cosmological term, together with other 
"modified gravity" terms, such as a Weyl curvature conformal gravity term. In 
addition, depending on the possible introduction of a fiber over spacetime given by 
a finite noncommutative geometry, one obtains further terms that give a Larangian 
for matter minimally coupled to gravity. This can be the Lagrangian of the minimal 
standard model (see [14]) or of the standard model with additional right handed 
neutrinos with Majorana mass terms (see [12j ) or models with supersymmetry (see 
[8] ) . In addition to the minimal coupling of matter to gravity one also finds terms 
such as a non-minimal conformal coupling of the curvature to the Higgs field. 

All this shows that one can use the spectral action as an action functional cither 
for pure gravity (on a commutative manifold) or for gravity coupled to matter 
on a product geometry. The main philosophy behind it is that pure gravity on a 
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noncommutative space can manifest itself as gravity coupled to matter when seen 
from a commutative point of view. 

The asymptotic expansion of the spectral action was also computed explicitly 
for truly noncommutative spaces like noncommutative tori, see (22j . [24j . 

We work here under the assumption that the spectral action is our modified grav- 
ity model, for commutative and noncommutative geometries alike and we discuss 
its behavior in a mixmaster case where the underlying spatial slices arc noncom- 
mutative 3-tori. 

Wc also discuss possible slow-roll inflation scenarios derived from the spectral 
action and how they affect and interfere with the underlying mixmaster dynamics. 

4.1. Slow-roll inflation in anisotropic cosmologies. In the usual isotropic 
Fricdmann cosmologies with Lorentzian metrics of the form 

ds^ = -dt^ + a{tf{dx^ + dy^ + dz^), 

inflation is an accelerated expansion of the universe that corresponds to the scale 
factor a{t) satisfying a > 0. The Friedmann equation relates the scale factor to the 
Hubble parameter 

a 
and the slow-roll models of inflation are based on the relation of the latter to a 
scalar field 4> with potential V{4)) via 

The slow- roll condition then corresponds to the condition that (fp' « |F(0)|, so 
that the term in (j) in the Friedmann equation becomes negligible. 

In the case of anisotropic spacetimcs of the form (12.11) , one introduces an average 
scale factor 

(4.1) a(i) ^ {a{t)b{t)cit)f'^ 

and directional Hubble parameters 

Hi = -, H2 = -, Hj, = -, 
a c 

and an average Hubble parameter 

H=^{Hi+H2 + Hi). 

This satisfies 

\ a b c 1 (dbc + bac + cab) 1 a 

" 3^a "^ 6 "^ c^ " 3 (a6c)2/3 (afec) i/3 = a' 

Thus, we obtain the same picture as in the isotropic case, but for the average 
Hubble parameter and the average scale factor. 

In the case of a mixmaster cosmology (|2.2p . where pi +P2 +P3 — 1, the average 
scale factor is just given by a = (tPi+P2+P3)i/3 ^ ^1/3^ ^-^^^ ^/^ ^ {l/3)t~\ Thus, 
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the Friedmann equation for the Hubble parameter in a pure mixmaster dynamics 
is of the form 

(4.2) ^- = H^\t-\ 

a 6 

or equivalently Hi + H2 + H^ = i~^. 

For a mixmaster cosmology, one also has 

a b c 

- + - + - = 0, 

a b c 

since, for a{t) = t^^, one has a/ a = pi{pi — 1)^^^, and similarly for the scale factors 
b{t) and c(i) so that 

- + 7 + - = (Pibi - 1) +P2iP2 - l)+P3iP3 - l))i~', 
a b c 

which vanishes, since we assume pi + P2 + Ps = ^ and pf + P'2 + P'i = ^■ 



4.2. The nonperturbative spectral action for 3-tori. Since we are dealing 
with spectral triples that are isospectral deformations of commutative tori, we can 
refer to the computation of the nonperturbative spectral action obtained in |31| for 
3-dimcnsional flat tori, using the Poisson summation formula as in |11) . 

Proceeding as in Theorem 8.1 of J3T], for a torus T^ with metric a{t)'^dx'^ + 
b{t)^dy^ + c{t)'^dz'^ and Dirac operator 



the spectral action is of the form 

A3 

'47r3 



(4.4) TT{f{fJA'^))^a{t)b{t)c{t)^^ I f{u^+v^+w^)dudvdw+0{A-''), 



for arbitrary fc > 0. 

Thus, for the mixmaster torus T^ with Dirac operator ^t as in (|2.9p . one finds, 
independently of the spin structure, 

(4.5) Tr(/(^,VA2)) =t--^ f f{u^ +v^+w^)du dv dw + 0(A-'=), 

47r3 7k3 

since in this case Pi+P2+P3 = 1, so that a{t)b{t)c{t) = ^pi+p^+ps ^ i. 

Notice that the factor a{t)b{t)c{t)K^ /{Att'^) in (|4.4p behaves exactly like the 
isotropic case, when one introduces the average scale factor a{t) — {a{t)b{t)c{t)Y^^ 
as discussed in the previous subsection. In terms of the average scale factor the 
spectral action has the form 

(4.6) Tr(/(^^/A')) = ^^^^"^f I f{u^+v^+w^)dudvdw +0{A-''). 

47r'^ ./b3 
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4.3. Noncommutative mixmaster cosmologies and inflation. We can then 
adapt the same analysis used in [31] , [32] , [9] to obtain a slow- roll potential out of a 
perturbation of the spectral action. To that purpose, we compute the spectral action 
in 4-dimensions, on a product T^ x S^, with a compactified direction S^ of size /3 
(after passing to a Wick rotation to Euclidean space and a compactification), with 
the product Euclidean metric. The replacement D^ M> D'^ + cj)^ produces a potential 
for the field 0, which, for sufficiently small values of the parameter x = (jP /A'^, 
recovers the usual shape of a quartic potential for the field (f>, conformally coupled 
to gravity. This gives, as in |5T| . 

(4.7) Tr{h{DyA')) = ^'/^°WH^)^(0 |°° ^/,(^) ^^ + o{A-'^y 
The perturbation is then computed as 

(4.8) TvihiiD^ + 02)/A2)) = TrihiD^/A^)) + ^^!M^M£W ^(^2/^2^^ 

An 

with V(a;) = /„ u{h{u + x) — h{u)) du, for x — (jp / A"^ where the last term 

A^PamtW) , , 
47r 
determines the slow-roll potential V{x). 

For a Kasner metric, this term is of the form 

with a{t) the average scale factor. This calls for a discussion of the i-dependcnce 
of the other parameters, A and /?, in this expression. 

4.4. Time dependence of the parameters. We have argued in [3T] and [35] 

that in the case of an isotropic expanding cosmology with a single scale factor 
a(i), the energy scale associated to the cosmological timeline should behave like 
A{t) ^ \/a{t). The possible interpretations of a time dependence for the parameter 
/? are less easily justified, as this parameter is an artifact of the S'^-compactification. 
We refer the reader to section 3.1 of [32] for a discussion of the interpretation of 
this parameter as an inverse temperature and its relation to the parameter A. 

In view of the expression of the spectral action (|4.6p . in terms of the average 
scale factor a(i) = (a(i)&(i)c(i))^/^, it seems natural to apply the same reasoning 
on the time dependence of the parameter A and /3, as in [32], in terms of a{t). This 
means, for example, that one expects the energy scale A, along the cosmological 
timeline, to behave like A(i) ~ l/a{t). 

In the usual setting of isotropic Friedmann cosmology, one relates the scale factor 
and the Hubble parameter to a density function, 

(4.10) (^)'-^^ = ^^W' 

where the form of the density function depends on the various cosmological eras: in 
a modern matter-dominated universe, the density function behaves, as a function 
of time, like p{t) ^ a{t)~^, which gives an evolution of the scaling factor with power 
law a{t) r^ t'^'^; in a radiation-dominated universe the density function behaves like 
p{t) ~ a{t)~^, and consequently the scaling factor evolves as a{t) ^ t^/^; while in 
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a vacuuni-domiiiated universe, as in the inflation epoch, p{t) is constant and a{t) 
is growing exponentiaUy. In a slow-roll model of inflation, the constant is given by 
the plateau value of the slow-roll potential. 

The latter observation implies that, in the expression (|4.9p for the potential 
obtained from the spectral action, assuming as above that A and a have an in- 
verse dependence on time, the consistency with (|4.10p would then suggest that the 
parameter /3, introduced artificially in the model as a radius of Euclidean compact- 
ification, should have a time dependence that is determined by the behavior of the 
density function /3(t) ^ p[t) (up to a multiplicative constant). 

As we will see in i i4.6l below, it is convenient to consider different possible t- 
dcpcndences of the parameter /3, as these make it possible to model transitions 
between different regimes in the very early universe (see also j30j for other discus- 
sions of time dependence of parameters in the very early universe in cosmological 
models based on the spectral action). 

Notice that, in any case, the time dependence (and the dependence on /3 and 
A) is only in the amplitude multiplicative factor of (|4.9|) . and it disappears entirely 
when one computes the slow-roll parameters, which depend only on the ratios V /V 
and V" /V (see [31] and [32]), through the expressions 

^ ~ 167r V V' / ' ^^ Sir V ' 



4.5. The scalar field potential and the mixmaster dynamics. In the presence 
of a slow-roll potential, an initial mixmaster dynamics gets altered: the coupling of 
the potential to the average scale factor via the (anisotropic) Friedmann equation 
disrupts the mixmaster oscillations. This phenomenon was already well known in 
the case of the classical mixmaster dynamics (see for instance chapter 8 of [35]). 
where it is used to argue that the potential provides a mechanism of transition from 
an anisotropic chaotic system to a standard big bang singularity. 

We discuss here the effect of the slow-roll potential associated to the spectral 
action on the underlying mixmaster dynamics, and on the corresponding transfor- 
mation of the torus moduli. 

It is easy to see why, already in the classical case, the presence of a slow-roll 
potential disrupts a mixmaster dynamics. In fact, suppose given, at time tg, an 
initial condition given by the Kasner metric of a Kasner epoch in a mixmaster 
universe. We have a(to) = t^^, b{to) — t^'^, c(to) = ^o^. Suppose that a slow-roll 
potential is turned on, with value near the plateaux level. Then the Friedmann 
equation predicts an evolution of the average scale factor by a{t) = o(fo) cxp(j(t — 
to)), where 7 = ^(87rGyoo)/3, with Voo the plateaux value of the potential. This 
is compatible, for example, with a solution for the individual scale factors of the 
form 

a{t)= iyV-i/3exp(7(i-io)), 
(4.11) b{t)^ 4/V-i/3cxp(7(f-fo)), 

c{t)^ iyV-i/3exp(7(f-fo)), 
where the exponential factor dominates over the underlying mixmaster dynamics. 
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4.6. Damping effects. It would be interesting to further study if one can use the 
presence of a slow-roll potential associated to the spectral action to model a tran- 
sition from a noncommutative early universe to a commutative modern universe, 
through a simultaneous "override" of the mixmaster oscillations both for the metric 
parameters and for the noncommutative torus moduli. 

In fact, according to the general philosophy about noncommutativity in the 
spacetime directions, one expects it to appear only close to the Planck scale (or 
the Planck era in the cosmological timeline), while near the unification scale (and 
therefore during the inflationary epoch, which is located between the unification 
and the electroweak epoch), the universe is already exhibiting noncommutativity 
only in the extra dimensions of the finite geometry that describe matter and forces, 
but not in the spacetime directions. Thus, one expects to find some mechanisms 
that damp the noncommutativity before the slow-roll inflation mechanism becomes 
relevant. 

For example, one can conceive of the following type of toy model scenario, where 
one obtains in two steps a transition from a mixmaster dynamics with noncom- 
mutative tori in the very early universe to a slow-roll inflationary cosmology with 
ordinary commutative tori. The first step is a deformation of noncommutative tori 
to commutative tori, obtained via a deformation of the Kasner exponents, which 
disrupts the mixmaster dynamics but leaves the Friedmann equation unaltered, and 
then a second step, in a purely commutative setting, where a deformation of the 
parameter j3 affects a second transition from a regime based on the Friedmann equa- 
tion for a Kasner metric to the Friedmann equation for an inflationary cosmology 
with slow-roll potential determined by the spectral action. 

4.7. Deforming the Kasner parameters. Suppose that the dependence of the 
exponents Pi{u) of ()2.3|1 on a parameter u is modified from a pure mixmaster 
dynamics, where u varies over the discrete set m„ determined by the corresponding 
dynamical system, to a continuous deformation with u > 0, with corresponding 
matrices 



l{u) = (^ l" J ) e Gum, 



replacing the sequence of matrices 7„ G GL2(Z) of (|2.1ip . The action of 7„ on 
the torus moduli, in either p.7|) or p.lO|) . then extends to a similarly defined 
action of the matrices 7(u) G GL2(IR), which is no longer by isomorphisms of ^e 
(respectively, Morita equivalences of Ae^^ ) but by transformations that deform ^e 
(respectively, Ae^) to a family of non-isomorphic Aq^u) (respectively, non- Morita 
equivalent Ag^(^u))- 

We can check the effect of such a deformation of the Kasner parameters in (j2.3p . 
with u = u{t), on the Friedmann equation. We find 

n —11 , ii. (u^ — 1) 



a l + u + vP- (1 

b 1 + M _-^ uu(u-\-2) 



t^' ,. . . , ,..;. log^ 



b 1+U + U^ (l+U + u2)2 

c u{l + u) _j u{2u + l) 

C ^ 1+U + m2 (l+y + u2)2 
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SO that the Hubble parameter 

d 1 ,a b c. 1 , 

a 6 a b c 6 
remains the same, along the deformation, as the Hubble parameter for the initial 
mixmaster dynamics (|4.2|) . In fact, the first terms in the above expression for the 
logarithmic derivatives of the scale factors add up to t~^, while the second terms 
add up to zero. 

Thus, through a deformation u = u(t) of this sort, which modifies the Kasner 
exponents, one can disrupt the mixmaster dynamics and deform the underlying 
noncommutative 3-torus ^e to a commutative torus, through the corresponding 
action on the torus moduli of the transformations 7(u) S GL2(IR) as above, while 
maintaining the Friedmann equation unaffected. 

4.8. Transition to slow-roll inflation. As a second step, after interrupting the 
mixmaster oscillation and damping the noncommutativity as described above, we 
want to obtain a transition from the Friedmann equation for a Kasner metric to 
the Friedmann equation for an inflationary cosmology. This can also be obtained 
through a deformation, this time of the parameter (3. 

In fact, for a fixed plateaux value Voo of the slow-roll potential V{x) of (|4.9p . we 
consider the expression 

Att ^°°- 

Here we may assume, as discussed in ^T^^ that A(t) ^ a{t)~^, which leaves us with 

mm,, mvoo 



An o(i) 47r ' 

A deformation from an initial phase with 

^ ' 3Voo t 
to a later phase with (3{t) ~ a{t) would then give the desired transition from 
the Hubble parameter of a Kasner metric to the one of a slow-roll inflationary 
cosmology. 



The two steps described in §473 and in this subsection arc only a toy model, 
since we do not provide a viable physical mechanism that produces the desired 
deformations affecting the transition to commutativity and to the slow-roll infla- 
tion scenario. We only exhibit a geometric model of how such transitions may be 
possible. 

5. Other aspects of the model 

In this final section we outline briefly other aspects of the model of mixmaster 
dynamics on noncommutative 3-tori, in relation to recently developed methods for 
explicit computations of the spectral action |10| : to some existing computations, in 
the commutative case, of the modified gravity terms in the asymptotic expansion of 
the spectral action for a Kasner metric [36] ; to the role of the diophantine condition 
in the asymptotic expansion of the spectral action for noncommutative tori [22], 
[24] : and to the effect of coupling gravity to matter [9]. We also mention other 
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physical models in which mixmaster dynamics play an important role and for which 
the model presented in this paper may have some relevance. 

5.1. The 4-dimensional Dirac operator for Euclidean Kasner metrics. We 

now look at the 4-dimensional Euclidean version of the Kasner metric, of the form 
dt^ + a{t)^dx^ + bitydy"^ + c{t)'^dz'^, and in particular at the case where a{t) = t^^ , 
b{t) = tP^, c{t) = tP^, with pi + p2 + p3 = 1 = pI + P2 + pI as above. We adapt 
some of the results obtained in [10] for the isotropic Euclidean Robertson- Walker 
metrics to this non-isotropic case. 



One can consider the Euclidean version of the anisotropic metric p.2p . of the 
form 

(5.1) dt'^ + a(tfdx^ + b{tfdy^ + c{tfdz^. 

One can proceed as in [lOj . and construct the corresponding 4-dimensional Dirac 
operator (with summation over repeated indices) 

(5.2) ip = Ye^^^^+\i^u.,ra'\ 

where 7*" are the gamma matrices, and the spin connection Usri now has nonzero 
terms 

a h c 

<^Wl = —, W202 — T: '^303 = -• 

a c 

Thus, one obtains an operator of the form 

with o the average scale factor as in (|4.ip . and where the operator D is given by 

,1 a .Id ol 9 

a ox oy c oz 

so that, as in the case of [TO], one has 7°!? = ^^t^ © ^t'*; since 7*'7-' = i(jj, with 
a sign difference due to our use of —icTj instead of iaj in (|2.6p . We also have, as in 
[TO], that 7°i:> = -D-i° so that D^ = ij°D)^. 

It is convenient here to rename the operator i^rps using the following notation 

(5-4) ^a,b.c ■■= -*(^l-^ + ^2T^ + '^3-^). 

(1 ux u (Jy c c/4. 

One can write the square of the Dirac operator ]p^ as 

(5.5) -d + ^-r+fa.b.c-^^^^, 

since the cross term contributes only a term of the form ^^M^D -\- D^^-S^ which 
gives 

, d .1. d 9,1,9 9,1.9. 

-*("^9^^a)9^ + "^9^^6)9^+"^9^(c)9^)' 

and, as above, 7°D = -^a,b,c ® ^a,b,c- 

Thus, one can reduce the spectral problem for this operator to an infinite family 
of one-dimensional problems, similarly to what happens in [TO]. Here one uses the 



(5.8) '{jr, + 7i-? + >?a,b.c->^e 



18 CHRISTOPHER ESTRADA AND MATILDE MARCOLLI 

fact that the two operators t^a,b,c and ^^ 52 ^ are simultaneously diagonalized 

with spectra 

(Kr\ G la \ r_Lo ii/^" + -Si/2 TO + S2/2 71 + 53/2 ^ ^ -^vSi 

(5.6) Spec(^a,6,c) = {±27r ( , r , ) {k,m,n) £ Z'^}, 

a c 

(5.7) 

Q (fH \ no ii Xfc+5i/2) b(m + 52/2) 0(^ + 53/2) ^,,,^,^ 3 

Spcc(^„2 ,2 ,2) = {±27r ( , , )\\\[k,m,n) gZ"}. 

"3"'-b"'~ a^ 0^ c^ 

We obtain in this way the one-dimensional problems 

d 3d 

where we write Xa,b.c for the elements in the spectrum (|5.6p and Xa2_ t^ c2. for those 

in the spectrum (|5.7p . This family (J5.8I) of one-dimensional problems is then, in 
principle, suitable for a computation of the spectral action based on the Fcynman- 
Kac formula, as described in [TOl . 

In the mixmaster case, one has spectra 

Al,, = t'''' (fc + y )' + 1'"'' (- + y )' + *-'^^' (" + y )' 

d ' b ' c ^ ^ ^ 

since a(i) = tv\ h[t) = tP\ c{t) = fP^ , with a/a'^ = pifP^-^r'^P^ = pii^^Pi+i) ^^^^ 
while the first term in the operators (j5.8p of the one-dimensional problems takes 
the form 



'^dt^ 2t' ' 



since o = t^/^ and d/a = t ^/3. 



5.2. Asymptotic expansion and the Kasner metrics. In the asymptotic ex- 
pansion of the spectral action over a commutative or almost-commutative geometry, 
one finds gravitational terms of the form (see [12] ) 

/ (i^ + yC^M-paC^"''" + ToR*R* - ^oR\hA dv, 

where Cf^^pa is the Weyl curvature, R* R* is the form representative of the Pontr- 
jagin class, which integrates to the Eulcr characteristic, and H is the Higgs field, 
with the coefficients expressed in terms of the momenta of the test function / in the 
spectral action and the Yukawa coupling matrix for the matter part of the model. 
An analysis of the gravitational terms in the asymptotic expansion of the spectral 
action in the case of a Kasner metric (and more generally for Bianchi V models) was 
given in [36] . In that paper, the authors focus on the modified gravity term given by 
the Weyl curvature that appears in the asymptotic expansion of the spectral action, 
neglecting the term with the conformal coupling of the Higgs field to gravity. They 
obtain modified Einstein equations of the form 

R""" ~ Ig'^R - o^o^l (2G'^,^:« - C'^'-'^R^^) = nlTZtter- 

They compute this explicitly in the case of Kasner metrics and other anisotropic 
Bianchi models, and they find that for such models the obtained equations indeed 
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differ from the ordinary Einstein-Hilbert ease, henee distinguishing the modified 
gravity model given by the spectral action from ordinary GR. 

In the mixmaster case, with ordinary commutative tori T^, one can see that 
their computation gives a correction term in the above equation with respect to the 
usual Einstein-Hilbert case of the form 

— iaoKn V-^ / // x2 N 

—^^4 2^Pr[PlP2P3 + Pi+l[(Pr - Pi+l) - PiPi+l) 

i 

+ iPi - ^)[-^Pi+liPi+l - 1) + ■^Pr+2{Pt+2 - 1) -PriPi - l)j 
+ ip't + 2 - 3k + (1 - p,){pi - ft+i - P,;+2) j (2pj - _p,;+i - K+2) j ■ 

As observed in |36j , this contribution is relevant in the early universe and becomes 
negligible for later times. 

5.3. The asymptotic expansion and the diophantine condition. In the case 
where the ordinary tori T^ in the mixmaster model are replaced by noncommutative 
tori Tq , dealing with the asymptotic expansion of the spectral action becomes more 
delicate. 

The perturbative expansion of the spectral action for noncommutative tori was 
computed recently in [22], [24]. In these calculations, one considers the spectral 
action 

(5.9) Tr(/(^yA2)), where ^^ = ^ + i, with A = A + eJ AJ-^ , 

where e is the commutation sign J^ = e^ J of the real structure involution J and 
the Dirac operator and A is a self-adjoint one form A = ^^ ai[^, bi], with a.i and bi 
in ^e- While, as we recalled in the previous sections, the Dirac operator ^ on Te is 
inherited via isospcctral deformation from a commutative torus, the operator ^^, 
twisted with the inner fluctuations A as above, depends on the noncommutative 
torus itself, and the results of [22], [24l show that this dependence manifests itself 
in an explicit dependence of the spectral action Tr(/(^^/A^)) on refined number- 
theoretic information on the modulus Q of the noncommutative torus. The compu- 
tation of the spectral action in [22] is through the asymptotic expansion for large 
A, in the form 

(5.10) TrififJA^)) = E ^ ^' /l^^l"' + /(O) C#. (0) + O(A-i), 

/c(Edini Spcc+ 

where the sum is over points in the strictly positive part of the dimension spectrum 
and the integration ^|^^|~*^ can be expressed as a residue of the zcta function of 
^A at the pole s = k. 

In particular, a diophantine condition on Q plays a crucial role in the results of 
[22], [H]. This "badly approximable" condition is formulated as follows. A vector 
V in K" is i5-diophantine, for some (5 > 0, if there exists a C > so that 

(5.11) \v-q-2TTk\>C\q\-^, V9eZ"\{0}, Vfc G Z. 

The set of diophantine vectors is the union over 6 of the sets of (5-diophantine 
vectors. A matrix Q is diophantine if there is a vector v G Z"" such that Qv is 
diophantine. Almost every matrix is diophantine (with respect to the Lebesgue 
measure) . 
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The positive dimension spectrum for an 7i-dimensional noncommutative torus Tq 
consists of the points {1, 2, . . . , n}, which are all simple poles, under the assumption 
that the matrix ^6 is diophantine. Under this same diophantine assumption, the 
residues at the points of the dimension spectrum are computed explicitly in [22]. 
The top n-th term is the same as for the unperturbed Dirac operator ^; the {n — k) 
terms with k odd vanish, and the {n — 2) term also vanishes. 

Thus, in the case of a 3-torus, the only nonvanishing term agrees with the unper- 
turbed case, so that the resulting perturbative spectral action under the diophantine 
condition does not differ in form from the ordinary torus case. By comparison, the 
case of a 4-dimcnsional torus exhibits a more interesting phenomenon whereby the 
perturbative spectral action in the case satisfying the diophantine condition has the 
form 

47r2 
Stt'UA' - —f{0)r{F^-^F^,) + 0{K-% 

where F^^'^ is while for the commutative case one would have I^a = ^, since the 
inner fluctuations would all be equivalent to zero in that case, so the expansion 
under the diophantine condition is different in this case. 

The transformation (j3.7p of the parameters 6_ in the mixmaster dynamics via 
the matrix (j2.1ip preserves the diophantine condition, hence in a phase of pure 
mixmaster dynamics, the asymptotic expansion of the spectral action of the non- 
commutative tori given in |22| . |24j remains valid. However, in the presence of a 
damping effect that destroys the mixmaster evolution and transitions to an isotropic 
and commutative torus geometry, the parameters 9_ are deformed to zero along a 
transformation that no longer preserves the diophantine condition. In fact, along 
such a deformation, the parameters 9_ will hit infinitely many values that do and 
do not satisfy the Diophantine condition and the properties of the asymptotic ex- 
pansion of the spectral action may vary accordingly in a seemingly chaotic manner, 
until the final stage of commutative isotropic tori is reached. Even though for 3- 
tori, unlike the case of 4-tori, the unperturbed commutative spectral action ends 
up looking the same, one should still investigate what happens in the intermediate 
stages of the evolution of the parameter 9_ in between an initial anisotropic mixmas- 
ter noncommutative case satisfying the diophantine condition and a final isotropic 
commutative case, especially at all the intermediate values of 9_ that do not satisfy 
the diophantine condition. 

5.4. Coupling to matter. In addition to pure gravity, described by the spectral 
action on the mixmaster tori (commutative or noncommutative), one can have a 
nontrivial coupling to matter, by taking the product of the spectral triple with a 
finite noncommutative geometry F, which specifies the matter content of the model, 
as in [H] and [M]. 

In terms of the asymptotic expansion of the spectral action, we can follow the 
computations of [53]. We now consider the example of a product geometry of a 
noncommutative torus Tq and a finite geometry of the form 

(5.12) Ap^M^iC), -Hp^MgiC), Dp^O, Jp = Jg, 

where Jq{L{T)) = L{T*), for T an element in the algebra and L{T) its representa- 
tion on the Hilbert space. 
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It is also convenient to write the spectral geometry for the noncommutative torus 
in the form TL = "Ht, the GNS representation for the tracial state T{a) = ao on Aq, 
so that, r being faithful, this is the Hilbert space completion of Aq in the inner 
product (a, b) — T{a*b). Then the real structure is Jg : Hr -^ 'Ht, Jo{h) = h* , with 
Jq = Jo and, for L{a) an element of Aq acting by left multiplication on "Ht, 

JoL{a)Jf^^h = JoL{a)Joh = J^L{a)h* = hL{a*) = R{a*)h, 

where R(a) is the action of a by right multiplication. Then the spectral geometry 
from the noncommutative torus can be written as 

A^A^, n = nr®C'^^"''\ ^^iSf,(g)-ft', J ^ Jo® Co, 

where C^ is an irreducible representation of C/^'*"^(M"), the 7'' are the gamma 

matrices implementing the representation, and Cq is the antiunitary operator on 
C^ implementing the charge conjugation. The (5^ are the basic derivations on 

the noncommutative torus, acting as 

^p I 2^ "'•"■ I - * 2^ rf^QrU- , 



(^a,uM =1^; 



on elements of T-Lt- written as a = "^^ arvT , where r G Z", and {a,.} G iS(Z"), and 
u^ is a Weyl element in Aq defined by 



exp 



j<k 



u, u^- ...U^ 



The sign of the spectral triple is 1 ® 7"+^ for n even. One has the relations 
JV = eVJ, and J^ = (' = (-1)L»/2J , The first implies that Co7^ = -e7^Co. 

As shown in |41j , this way of writing the spectral triple structure on the noncom- 
mutative torus is equivalent to the one we used above by isospectral deformation. 

The product geometry then has algebra AQ<E)Mq{C), Hilbert space H® Mq{C) — 
Ht <E) C^ (8) Mg(C), Dirac operator D = i5^ 7'' idg and real structure 

J = Jo(E)Co(g)Jg. 

Then, following the same computation of [24] , one finds that the presence of the 
finite spectral triple F of (|5.12p alters the perturbative expansion of the spectral 
action by a factor of q^ in the leading order term. For example, in the case of a 
4-dimensional torus, under the assumption that the diophantine condition on Q 
holds, one finds 

47r2 

In the case of the commutative mixmaster tori T^, one can also proceed as in [S] 
to compute the effect on the gravitational part of the spectral action of the matter 
sector, which as in [9] also delivers an overall multiplicative factor q^, equal to the 
rank of the finite geometry. 

Thus, in the case of the mixmaster tori, switching back to the same notation 
used in the previous sections, we find that the spectral action, for this case with 
the finite geometry F of (j5.12p . 

q^aitfA^ f ^,^,2 , 2 , 2N 



(5.13) TT{f{D^/A'))^^^-^ / fiu' + v'+w')dudvdw. 
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Correspondingly, when one computes the inflation potential as in (|4.8p . one finds 
that, as in the case of [9], the resulting slow roll potential is q'^V{x), with V{x) the 
inflation potential in the absence of the finite geometry. Since only the amplitude 
of the potential is affected, the slow-roll parameters do not change, but as argued 
in [32] and [9], the amplitudes for the power spectra for density perturbations and 
gravitational waves (scalar and tensor perturbations) detect the different scaling 
factors in the slow-roll potentials and can therefore, in principle, detect the presence 
of the matter sector through the rank of the finite geometry. 

5.5. Relation to other physical models. We mention here briefly other recent 
cosmological models where the kind of noncommutativc deformation of the mix- 
master dynamics we described above may turn out to be useful. 

A first possible context is Hofava-Lifschitz gravity. It was recently shown in 
[3] that the mixmaster universe provides a mini-superspace truncation of the field 
equations of Hofava-Lifschitz gravity. The latter is a recently introduced higher 
derivatives modified gravity theory, applicable in the ultraviolet regime |25| . The 
Lagrangian density for this theory is derived from a superpotential that contains a 
Chern-Simons gravitational term and a 3-dimensional Einstein-Hilbcrt term. The 
resulting expression has terms of the form aR+fS+jCijC^^ +SCijR^^ +eRijR'^^ +(R^, 
where the first two terms give rise to the usual Einstein-Hilbcrt action and the 
remaining higher derivative terms contain the Cotton tensor C^^ (see (2.16) of [5]), 
which vanishes if the 3-dimensional spatial sections are conformally flat, and the 
curvature tensors. The action is invariant under coordinate transformations of the 
spatial sections. 

Another such context is loop quantum cosmology. It was recently shown in [40] 
that, in the setting of loop quantum cosmology, one finds an oscillatory behavior of 
mixmaster type as one approaches the singularity as a simplified system from the 
equations of motion in the Hamiltonian formulation (sec §2.3-2.8 in [40]) and that, 
moreover, on the resulting reduced phase space one can also introduce a scalar field 
with an inflation slow-roll potential (see §2.7 and chapter 4 of [10]). The presence 
of the scalar field has a damping effect on the mixmaster oscillations in this model 
(see Appendix A of [40]). 
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